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Graphene is a two dimensional crystal of carbon atoms with fascinating electronic and morpholog- 
ical properties. The low energy excitations of the neutral, clean system are described by a massless 
Dirac Hamiltonian in (2+1) dimensions which also captures the main electronic and transport prop- 
erties. A renormalization group analysis sheds light on the success of the free model: due to the 
special form of the Fermi surface which reduces to two single points in momentum space, short range 
interactions are irrelevant and only gauge interactions like long range Coulomb or effective disorder 
can play a role in the low energy physics. We review these features and discuss briefly other aspects 
related to disorder and to the bilayer material along the same lines. 

PACS numbers: 



I. INTRODUCTION. 

For the last few decades, the Landau Fermi liquid (LFL) has been the standard model in condensed matter to 
describe most metals [TJ [2J . A renormalization group (RG) analysis of the effective continuum model of non-relativistic 
electrons allows to understand this universal behavior as a consequence of the marginal character of electron-electron 
interactions in the presence of a Fermi surface [66] . Strongly correlated materials such as heavy fermions or the hight 
temperature superconductors have given rise to very important new ideas in the last 20 years due to their "non- 
Fermi liquid" behavior that arise from their strongly interacting nature [3|. The appearance of graphene in 2004 [3] 
has originated a "second revolution" in the field due to its very many exotic electronic and morphological properties. 
Unlike the cupratcs where there is not up to today an agreement on the model to describe the basic material, graphene 
resembles the situation in quantum field theory (QFT) where the starting point is clear: the low energy excitations 
can be modeled by the massless Dirac Hamiltonian in two spatial dimensions. Under the RG point of view, graphene 
has defied the LFL paradigm by confronting us with a material described by a universal non-interacting Hamiltonian 
that shows agreements and contradictions with the usual LFL behavior. From many other points of view it is a 
fascinating system which lives in between different branches of physics as condensed matter, QFT, statistical physics 
or the theory of elasticity. 

Graphene is a two-dimensional crystal of carbon atoms arranged in a honeycomb lattice: a single layer of graphite. 
Its synthesis 016], amazing properties (it is flexible like plastic but stronger than diamond, and it conducts electricity 
like a metal but is transparent like glass |67j ) and potential applications [5], have granted the 2010 Nobel prize in 
Physics to their fathers A. Geim and K. Novoselov. 

One of the most interesting aspects of the graphene physics from a theoretical point of view is the deep and fruitful 
relation that it has with quantum electrodynamics (QED) and other quantum field theory ideas [5HT2"]. The connection 
arises from the mentioned fact that the low energy excitations can be modeled by the massless Dirac equation in two 
spatial dimensions a fact that makes useless most of the phenomenological expressions for transport properties in 
Fermi liquids. From the morphological point of view, the mere existence of two dimensional crystals has been argued 
not to be possible [7j . Moreover despite of being one of the most rigid materials with a Young modulus of the order 
of terapascals, the samples show ripples of various sizes whose origin is unknown. The intrinsic relation between the 
morphology (honeycomb lattice) and the electronic properties constitutes a great playground for physical ideas as 
well as one of the sources of potential applications as a single molecule detector. 

Under the QFT point graphene has given rise to very interesting developments: the so-called axial anomaly [§1 HO] 
has acquired special relevance in relation with the recently discussed topological insulators [131 114) which provide a 
condensed matter realization of the axion electrodynamics [T3 HI] • Charge fractionalization has also been explored 
in the honeycomb lattice with special defects [17j [18] and quantum field theory in curved space [19l [20] and cosmo- 
logical models [5TJ [23] have been used to explore the electronic properties of the curved material. A very interesting 
development is associated to the generation of various types of vector fields coming from the elastic properties or from 
disorder that couple to the electrons in the form of gauge fields [23] . 

Among the unexpected properties related with interactions on the condensed matter point of view are the anomalous 
behavior of the quasiparticles decaying linearly with frequency |24j , and the renormalization of the Fermi velocity at 
low energies [35] . These properties were predicted theoretically in the early works [TTJ [2S] from the RG analysis that 
will be described in this work. We will review the RG aspects of graphene in the presence of Coulomb interactions with 
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special emphasis in the similarities and discrepancies with QED(3). We will show that the system has a non-trivial 
infrared fixed point where Lorentz covariance is restored appearing as an emergent property [27 . The RG analysis 
allows to classify graphene as a new type of electron liquid described by a Lorentz covariant infrared fixed point whose 
effective coupling constant is the fine structure constant. We will focus on the general, conceptual features referring 
to the specific bibliography for the technical details. The differences and similarities with the usual QFT models will 
also be emphasized. 



II. CONTINUUM MODEL FOR THE LOW ENERGY EXCITATIONS OF GRAPHENE 



The graphene structure is made of carbon atoms located in a honeycomb lattice. Three of the four available 
electrons of the carbon atom: 2s, 2p x and 2p y orbitals hybridize in a so-called sp 2 , a strong covalent bond that 
"draw" the honeycomb lattice and have typical energies around 3 eV. The fourth 2p z orbital perpendicular to the 
plane formed by the sigma bonds remains delocalized and is responsible for the metallic properties of the material. 
The a bonds give rigidity to the structure, while the 7r bonds give rise to the valence and conduction bands [55j. The 
electronic properties around the Fermi energy of a graphene sheet can be described by a tight binding model with 
only one orbital per atom, the so-called 7r-electron approximation and were obtained in the early works 29, 30J. The 
nearest-neighbor tight binding approach reduces the problem to the diagonalization of the one-electron Hamiltonian 



n = E a t a J 



where the sum is over pairs of nearest neighbors atoms i, j on the lattice and a%, are the usual creation and 
annihilation operators. What is special about the honeycomb lattice shown in Fig. [I] is the fact that, even though all 
carbon atoms are equal, their positions in the lattice makes them topologically non-equivalent. It is easy to see that 
to generate all the points in the lattice we must take a pair of (red and black) points. The lattice can be seen as two 
interpenetrated triangular sublattices of red and black points or, more technically, it is a lattice with a basis. This 
special topology of the honeycomb lattice is at the heart of all the electronic peculiarities of graphene. The Bloch 
trial wave function has to be build as a superposition of the atomic orbitals from the two atoms forming the primitive 
cell black (A) and red (B): 

* k (r) = C A 4> A + C B <f> B - (2) 
The eigenf unctions and eigenvalues of the Hamiltonian are obtained from the equation 
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= E(k)[ ri ) - ( 3 ) 
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where Uj is a triad of vectors connecting an A atom with its B nearest neighbors and vj the triad of their respective 
opposites, a is the distance between carbon atoms and e is the 2p z energy level, taken as the origin of the energy. The 
tight binding parameter t estimated to be of the order of 3eV in graphene sets the bandwidth (6eV) and is a measure 
of the kinetic energy of the electrons. The eigenfunctions arc determined by the the coefficients Ca and Cb solutions 



FIG. 2: (Color online) The dispersion relation of graphene. The six Fermi points of the neutral material are located at the 
vertices of the hexagonal Brillouin zone. 



of equation pL The eigenvalues of the equation give the energy levels whose dispersion relation is 




^3 \^3 3 

-E'(k) = ±t\l 1 + 4 cos 2 —^-ak x + 4 cos —^-ak x cos -ak y , (4) 

in which the two signs define two energy bands: the lower half called the bonding 7r band and the upper half called the 
antibonding n* bands, which are degenerate at the K points of the Brillouin zone. The dispersion relation is shown 
in Figure [2] Within the it electron approximation each site of the graphite honeycomb lattice yields one electron to 
the Fermi sea and the band is at half-filling. Since each level of the band may accommodate two states due to the 
spin degeneracy, and the Fermi level turns out to be at the midpoint of the band, instead a whole Fermi line, the 2D 
honeycomb lattice has six isolated Fermi points which are the six vertices of the hexagonal Brillouin zone. Only two 
of them are inequivalent and can be chosen as K\ — 0) and Ki = —K\ 31 . 

A continuum model can be defined for the low energy excitations around any of the Fermi points, say K\, by 
expanding the dispersion relation around it: k = K\ + (5k what gives from ([3]) the Hamiltonian 

n--\ta(, h ° ^ + ^V (5) 



2 I 8k x — iSky 

The limit lim a _>.o H /a defines the continuum Hamiltonian 

H = VF^.k, (6) 

where a are the Pauli matrices and the parameter vp is the Fermi velocity of the electrons estimated to be i>f ~ 
Ma ~ c/300. Hence the low energy excitations of the system are massless, charged spinors in two spatial dimensions 
moving at a speed vp. We must notice that the physical spin of the electrons have been neglected in the analysis, 
the spinorial nature of the wave function has its origin in the sublattice degrees of freedom and is called pseudospin 
in the graphene literature. In the absence of interaction mixing the two Fermi points it can be shown that they are 
topologically stable [3 2) and can not be lifted by smooth continuum perturbations of the Hamiltonian, another way 
of seeing the RG results that will be analyzed here. The same expansion around the other Fermi point gives rise to a 
time reversed Hamiltonian: T-ii = vp(—a x k x + a y k y ). The degeneracy associated to the Fermi points (valleys in the 
semiconductor's language) is taken as a flavor. Together with the real spin the total degeneracy of the system is 4. 
The calculation of the density of states for the linear dispersion relation E(k) — Vp\k\ gives 

(x? \ 9S \E\ . . 



where gs =4 is the degeneracy. Eq. Q shows one of the peculiarities of graphene as opposed to a usual two 
dimensional electron gas (2DEG): the density of states grows linearly with the energy and vanishes at the Fermi energy 
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(in the standard 2DEG the DOS is a constant in two dimensions). This fact has very important phenomenological 
consequences that have been tested experimentally. For our purposes it implies that the Coulomb interactions between 
the electrons is not screened in graphene. 

We finish this section by mentioning that the massless spinorial nature of the graphene quasiparticles, the linearity 
of the DOS, the degeneracy of 4 and the estimated value of the Fermi velocity have been experimentally tested [7] what 
makes the present model universally accepted in the community. We also mention that many electronic and transport 
properties can be described by the non-interacting model and that disorder which usually leads to localization in two 
dimensions, does not seem to play an important role in the graphene system. 

III. RG ANALYSIS: TREE LEVEL 

As discussed in the article by R. Shankar in this volume, the main issue of the renormalization group ideas as applied 
to condensed matter systems is that for special systems (critical, renormalizable) the low energy physics is governed by 
an effective Hamiltonian made of a few marginal interactions that can be obtained from the microscopic high-energy 
Hamiltonian in a well prescribed manner [331 134] . Following the nomenclature of critical systems, interactions are 
classified as relevant, irrelevant or marginal according to their scale dimensions. These dimensions determine whether 
they grow, decrease, or acquire at most logarithmic corrections at low energies. The effective coupling constants of 
a model at intermediate energies can be obtained by "integrating out" high energy modes even if there is no stable 
fixed point at the end of the RG flow. The Luttinger and Fermi liquids arc identified as infrared fixed points of the 
RG applied to an interacting metallic system in one or more dimensions respectively. 

In this context, the universal Landau Fermi liquid behavior of most electronic systems can be established in the 
form of almost a theorem: a system of electrons will behave as a Fermi liquid if: (1) The spatial dimension is D > 2. 
(2) The system has a smooth extended Fermi surface. (3) The interactions are repulsive, non-singular and short 
ranged. 

The key ingredient in the result is the existence of a finite Fermi surface which 

• Ensures a finite density of states at the Fermi energy and hence the screening of the electron-electron interaction 
making it effectively short ranged. 

• Makes the tree level scaling analysis effectively one dimensional as only the momentum component perpendicular 
to the Fermi surface scales like the energy. 

• Provides the kinematics necessary for the Landau channels to work. 

Most of the non- Fermi liquid behaviors searched for in two dimensional systems in the last years have been associated 
to either specific shapes of the Fermi surface as in the so-called Van Hove scenario 35 or to singular interactions 
|36| . The case of graphene is specially interesting for the reduction of the Fermi line to two single points. This fact 
prompted the early RG analysis [IT] which at the time aimed to find non- Fermi liquid behavior coming from the fact 
that Fermi points were characteristics of one dimensional systems with different universal behavior (Luttinger liquid). 

Next we will classify the local interactions in the graphene system at tree level system following [33J IM] ■ From the 
Hamiltonian ^ we can construct the action 

So = J d 2 k duo * [(wl - vpcr.k] (8) 

what allows to fix the scale dimension of the fields to [&] = —2 in units of energy. Notice that for usual non 
relativistic fermions the dispersion relation is also linearized around the Fermi surface. The crucial difference between 
the two cases is that, because the Fermi surface is a point, there is an isotropic scaling of the momentum so that 
under a rescaling of the energy u> ~ su, d k ~ s 2 k. In the presence of an extended Fermi line, we would have 
d 2 k = dk\\ dk± ~ dk\\ sdk± and the dimension of the fields would be [$f] — —3/2 in any number of dimensions. 

Having established the different scale dimension of the fields it is obvious that four (or more) Fermi interactions - i. 
e. all local interactions in the condensed matter approach - are irrelevant. (Notice that they were already generically 
irrelevant also in the non-relativistic case and only the Landau channels corresponding to special kinematics around 
the Fermi line were made marginal). Since we can not find any marginal local interactions in the graphene case this 
finishes the analysis. No magnetic or superconducting instabilities driven from local interactions are to be expected 
in the clean graphene system assumed that the couplings are small and perturbation theory applies, a fact that also 
seems to be confirmed by most experiments. 

Although the only possible local interactions in non-relativistic electrons are made of polynomials in the fields we 
know that relativistic fermions interact through gauge fields and these interactions will be marginal. We will describe 
how they arise in the graphene system in the next section. 
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IV. COULOMB INTERACTION. GRAPHENE VERSUS QED 

The unscreened Coulomb interaction is usually written as 

n mt = f /A, fft> *M^l!W , (9, 

2tt7 7 |ri-r 2 | 

The non-local interaction ^ can be understood as an effective interaction coming from a standard gauge field coupled 
to the fermions through the minimal coupling prescription giving rise to the interacting action: 



S int =9 J d z x dt J»(x,t)A^(x,t) , (10) 

where g = e 2 /4ttvf is the dimensionless coupling constant, the electronic current is defined as 

= (VPT^, vf^ct^) • (11) 

Notice that although the electrons are confined to live in the two-dimensional surface, the Coulomb interaction among 
them lives in the three dimensional space. We then face the problem of coupling the two-dimensional current to a three 
dimensional gauge potential. In ref. [TT] this difficulty was solved by integrating the usual photon propagator over 
the z component what produces an effective two-dimensional propagator for the gauge field with a l/|k| dependence. 



Since by construction the vector potential scales like the derivative, the interaction ( 10 1 is marginal. The full action 



S = g / d 2 xdt^ [d^ + igA^ * , A* = 0,1, 2 (12) 



is scale invariant. 7 M are a set Dirac matrices constructed with the Pauli matrices. The action (12) looks like that 
of (non-relativistic) quantum electrodynamics in two spatial dimensions but there are important differences that 
affect the discussion of the one loop renormalization of the coupling. It is the anomalous photon propagator what 
makes the interacting graphene system different from QED(2+1) and the reason why the results found there can 
not be translated directly to the physics of graphene. The coupling constant of QED (2+1) has dimension of mass 



(VM) and the theory is called "superrenormalizable" . It means that it has less divergences than its four dimensional 



counterpart. In fact, there are no ultraviolet infinities in QED(2+1). By construction, the model given by (12) is fully 
scale invariant, the coupling constant is dimensionless and, in this sense, the system resembles more what happens in 
QED(3+1) which is strictly renormalizable - or at a critical point- . This is also what induces a renormalization of 
the electron self-energy in the graphene system not present in planar QED. 



V. RG AT THE ONE LOOP LEVEL 



Once we have identified the marginal character of the Coulomb - and in general any generalized gauge-type 



interaction of the electronic current (11) with a vector field of weight 1 - at tree level we have to establish the 
behavior of the renormalized coupling constant. 

From the structure of the perturbative series we see that the effective -dimensionless- coupling constant is 

e 2 
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Airvp 



where e is the electron charge and Vp is the Fermi velocity. This is the equivalent of the graphene fine structure 
constant where vf replaces the speed of light c. We will proceed as in the Fermi liquid analysis and perform a 
perturbative renormalization at the one loop level. This amounts to assume for the time being that the effective 
coupling constant is small, something which is not guaranteed in the suspended samples. We will comment on that 
later. (The bare value of this constant at energies around leV is of the order of 2-3, corresponding to 300aQ ED ). 



The RG analysis of the model ( 12 ) was performed in full detail in refs. [TTJ 123 [55] • We will here summarize the 
results obtained there. 

Due to the very special structure of the model discussed in the previous section, the coupling constant renormal- 
ization of the model at one loop level is more subtle than both cases, the non-relativistic condensed mater fermions, 
and QED. The building blocks of the Feynman graphs are the electron and photon propagators: 



(13) 



G 





FIG. 3: Diagrams renormalyzing the model at one loop 



n ° (k) = 57OT' (14) 

and the tree level vertex T M = ie(j°, v 7). 

Due to the non-covariant form of the electron propagator (131 the system has four free para meters, the electric 
charge e, the Fermi velocity v, and the electron and photon wave function renormalization. A Ward identity associated 
to charge conservation relates the electron propagator and the vertex function so the renormalization of the model 
can be done as in QED by renormalizing the electron and photon propagators only. 

The one loop diagrams are shown in Fig. [3| It can be shown to all orders in perturbation theory that the photon 
propagator is finite and hence the electric charge is not renormalized. All the renormalization of the system comes 
from the electron propagator. The renormalization of the Fermi velocity and the electron wave function are obtained 
from the electron self-energy X computed from Fig. [3] (a) and from the relations: 

E(w > k) = ^(«,k)[w7 -Z v (w,k)«7.k] . (15) 



The wave function renormalization 



defines the anomalous dimension of the field 



v = Z v v R . (16) 



~ — t^j — -|w=o 1 (17) 



j = dlog Zy/dl, (18) 
(I is the RG parameter) and, hence the asymptotic behavior of the fermion propagator: 



The computation gives 



G(u,k) -J-j . (19) 



Z„ = l-|logA, (20) 



from where we can relate the Fermi velocity at two different energies 

v(E) = v(E ) 



(21) 



The last equation means that the Fermi velocity grows as the energy decreases, the main result of the analysis. The 
anomalous dimension is found to be 



*=1^ < 22 > 

Summarizing, the results of |11| are the following: 

• The Fermi velocity grows as the energy is decreased. This result implies a breakdown of the relation between 
the energy and momentum scaling, a signature of a quantum critical point. 
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• The electron-photon vertex and the photon propagator are not renormalized at the one loop level. This means 
that the electric charge is not renormalized, a result that could be predicted by gauge invariance, and it also 
implies that the effective coupling constant g = e 2 /Airvp decreases at low energies defining an infrared free fixed 
point of the RG. 

• The beta function of the Fermi velocity has a non-trivial zero at the value v = c. Hence the infrared fixed point 
defines a weak coupling model ruled by the fine structure constant of QED which is Lorentz covariant. 

• There is a critical exponent which determines the universality class of the model. The anomalous exponent 
means that the fixed point is not a Fermi liquid and resembles more a Luttinger liquid. 



VI. MARGINAL FERMI LIQUID? 



Because of the smallness of the ratio vp/c and for simplicity the Coulomb interactions in graphene are often 
modeled with a static scalar field that couples only the the charge density. This model was analyzed in 37J. The 
original motivation of this work was to explore if the gapless free infrared fixed point encountered in within the 
perturbative RG scheme would survive in the simplest non-perturbative resummation consisting in computing the 
one loop electron self-energy with an effective interaction given by the bubble sum in the photon propagator (1/N 
approximation). The logarithmic renormalization of the electron wave function found in this case has led to some 
confusion on the marginal Fermi liquid nature of the system that it is worth to clarify. 

The so-called marginal Fermi liquid (MFL) behavior was set as a phcnomenological model to explain some 
anomalous behavior of the cuprates high-T c superconductors. It has two prominent features: a linear behavior of the 
scattering rate with the energy near the Fermi surface 

X(k,Lj)~ g 2 N 2 (0)[cj l °^-i~x}, (23) 

where x=max[|w|,T] , and, more important, a vanishing quasiparticle residue at the Fermi energy: 

G(k, u) = 1 = p— — + G incoh , (24) 



Z u =1 ^ L=E, -log — 



^ = l--^U^-log|^|. (25) 

The two are related by causality since they are extracted from the real and imaginary part of the electron Green's 
function. The quasiparticle weight Z% vanishes logarithmically at the Fermi surface Ek = and the Green's function 
is entirely incoherent. As pointed out in the original MFL reference, such a behavior is only possible if perturbation 
theory, starting from a non-interacting Fermi gas, breaks down at some energy. We must also notice that the 
postulated MFL behavior in the high- Tq superconductors is expected to arise from the strongly correlated nature of 
the electron-electron interaction in these compounds. In the case of graphene the linear scattering rate comes from 



the singular nature of the unscreened Coulomb interaction (iV(0) = in eq. (23)) 



The (instantaneous) Coulomb interaction in graphene is described by the scalar component of the gauge field: 

H int = e J d 2 r*(r)*(r)^(r) , (26) 

where the scalar field has the propagator 

0(ri,*i;r a ,t 2 ) = ^-5(ti-t 2 ) 1 1 , . (27) 

47T |ri — r 2 | 

Inserting the RPA sum in the electron propagator we get the divergent contributions to the electron self-energy from 
which we can extract the following RG equations [BS] 

dZy 8 ,„ 2 — q 2 arc cos q. 8 1 

r(2+ y - 7=4;) + —, (28) 



<9l0gA 7T 2 g y/l - g2 7T g' 

dv F 8 arc cos g 4 v F 

j v f( 1 + — u J + Z — > I 29 ) 



dl0gA 7T 2 gy /l - g 2 7T g 
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where A is the ultraviolet cutoff and g is proportional to the effective coupling constant g = ^ . The simultaneous 
resolution of these equations shows that the quasiparticle residue goes asymptotically to a non-zero constant value 

EH- 

The infrared behavior of the system is that of a strange Fermi liquid: although the scattering rate grows linearly 
with the energy, the wave function renormalization Z runs to a finite value at the Fermi point. This paradoxical 
behavior can be seen as another manifestation of the "incompleteness" of the static model which can not be used to 
establish the infrared nature of the system. 

VII. OTHER RG ASPECTS 

A. Disorder and interactions 

In graphene many classes of lattice defects can be described by gauge fields coupled to the two dimensional Dirac 
equation [53] . The standard techniques of disordered electrons [4jl] can be applied to rippled graphene by averaging 
over the random effective gauge fields induced by curvature or elastic deformations. A random distribution of defects 
leads to a random gauge field, with variance related to the type of defect and its concentration. There is an extensive 
literature on the problem, as the model is also relevant to Fractional Quantum Hall states [41] and to disorder in 
d-wave superconductors. A random field, when treated perturbatively, leads to a renormalization of the Fermi velocity 
which makes it to decrease at lower energies opposing the upward renormalization induced by the long range Coulomb 
interaction. The simultaneous presence of interaction and disorder gives rise to new interesting fixed points. The 
issue was analyzed in [42] . The most interesting case arises when considering a random gauge potential which models 
elastic distortions and some topological defects. There is a line of fixed points with Luttinger-like behavior for each 
disorder correlation strength A given by v F — 2e 2 /A. An extensive analysis of the issue in the large N limit is done 
in [43]. 

B. Short range interactions. 

As we have seen short range interactions, such an onsite Hubbard term U are irrelevant, a fact due to the vanishing 
density of states at the Fermi level. As mentioned above, the density of states at low energies is increased by the 
presence of disorder. This, in turn, enhances the effect of short range interactions that were discussed in the early 
work |38j and have attracted some attention recently. These interactions can be relevant in the strong coupling regime 
of the hexagonal electronic and optical lattices [44H46] . 

A summary of the situation following [35J is as folllows: In the absence of disorder, an onsite Hubbard term favors 
antiferromagnetism. An antiferromagnetic phase, however, is likely to be suppressed by disorder, especially by the 
presence of odd numbered rings in the lattice. Then, the next leading instability that such an interaction can induce 
is towards a ferromagnetic phase. 

If a magnetic phase does not appear, electron electron interactions, even when they are repulsive, will lead to 
an anisotropic ground state. The existence of two inequivalent Fermi points in the Brillouin zone suggests that 
the superconducting order parameter induced by a repulsive interaction will have opposite sign at each point. The 
corresponding symmetry is p-wave. Disorder, in addition to the enhancement of the density of states mentioned 
already, will lead to pair breaking effects in an anisotropic superconducting phase. 

C. Bilayer graphene 

The synthesis characterization and analysis of the bilayer material (BG) occurred simultaneously with its monolayer 
counterpart |47j . The band structure was also theoretically described in the early publication [30] . The most common 
structure is made of two graphene layers arranged in the so-called Bernal stacking as shown in Fig. [4] It was soon 
realized that the BG was even more promising than the monolayer due to the possibility of opening and controlling 
a gap in the system [35] [35]. Under a QFT point of view it constitutes an interesting example of a chiral system with 
quadratic dispersion relation. 

A tight binding approach done to the minimal Bernal stacking bilayer shows that the system has four bands: Two 
low energy bands touching at a Fermi point and two at higher energy. An effective low energy continuum model [50] 
keeping only the lower bands can be described by the Hamiltonian 
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FIG. 4: Bilayer graphene lattice. 



=~ ( ° k ; ) , (so) 

where k — (k x +ik y ). The dispersion relation is quadratic but still has a non-trivial chirality. Unlike the free monolayer 
system whose effective description corresponds to regular massless fermions, there is no QFT with quadratic propagator 
and Dirac structure. In this sense, although being chiral, the Fermi points of the bilayer can not be termed as Dirac 
points. The interacting model has been subjected to a great interest recently due to the enhanced quality of the 
samples [5TI 152] , 

The role of interactions and possible instabilities in the bilayer system was first analyzed in |53j and has received 
renewed attention more recently [54TI57] . RG studies have been performed in [581460] . The interest of this case relies 
on the fact that while the quadratic dispersion makes short range interactions marginal, the existence of Fermi points 
(although with finite density of states) can give rise to departures from the Fermi liquid behavior. Marginal Fermi 
liquid behavior similar to the one discussed for the monolayer in Sect, 
complete list of magnetic and pairing instabilities has been described in 



VI has been advocated in [58] and a fairly 



VIII. AN ATTEMPT OF A SUMMARY 



The infrared properties of graphene described in this work explain the success of the non-interacting model to 
describe most of the low energy electronic and transport properties of the system. We have seen that the only 
interactions that can play a role in the system are of the form of gauge couplings. Of these, the most important 
is the unscreened Coulomb interaction that runs in the infrared to a very weakly interacting fixed point. In the 
non-relativistic model the interaction is marginally irrelevant. 

As a very interesting theoretical issue that makes the model differ from both its QFT and its condensed matter 
counterparts is the fact that the coupling constant renormalization at the one loop level is due to the upward renor- 
malization to the Fermi velocity to the infrared. This very fact sets a lower bound on the validity of the widely used 
instantaneous model for the Coulomb interactions which ceases to be valid when the Fermi velocity approaches the 
speed of light. Although this bound is so small that can be ignored for most practical matters [52] it is conceptually 
important that the infrared nature of the system has to be decided with the full retarded model. 

Another interesting issue is that the running of the Fermi velocity stops at a non-trivial fixed point where its value 
equals the speed of light so that the Lorentz covariance appears as an emergent property [27j . Unlike most infrared 
stable QFT where the fixed point is trivial (g=0), the graphene system has a non-trivial fixed point characterized by 
the fine structure constant. 

Within the static approximation a two loops calculation of the electron inverse lifetime shows that it is linear in 
the energy, a behavior that persists in the 1 /N approximation and that is anomalous for a metallic system. 

Experimental indications of both the behavior of the inverse lifetime of the electron linear in the energy [2H 2331 and 
of the Fermi velocity renormalization [2 5) have been recently reported. The irrelevance of the short range interactions 
excludes a priori any intrinsic superconductivity or magnetism in the samples. The RG results also exclude the 
possibility of opening a gap in the weak coupling regime, a fatal situation for the electronic applications. The 
possibility of opening a gap in the spectrum with or without the help of a magnetic field is crucial for the possible 
electronic applications and was the main concern in the strong coupling analysis. A recent update of the issue and 
the main references can be found in [64 . A non-perturbative RG analysis of the system in the lattice has been done 
in [55] . 
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Notice that there is an erratum in equation (10) of ref. [37] where the derivative of \ogZ appears instead of Z. 



